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Abstract 

We derive the effective equations for the out of equihbrium time evolution of 
the order parameter and the fluctuations of a scalar field theory in spatially 
flat FRW cosmologies. The calculation is performed both to one-loop and in a 
non-perturbative, self-consistent Hartree approximation. The method consists 
of evolving an initial functional thermal density matrix in time and is suitable 
for studying phase transitions out of equilibrium. The renormalization aspects 
are studied in detail and we find that the counterterms depend on the initial 
state. We investigate the high temperature expansion and show that it breaks 
down at long times. We also obtain the time evolution of the initial Boltzmann 
distribution functions, and argue that to one- loop order or in the Hartree 
approximation, the time evolved state is a "squeezed" state. We illustrate 
the departure from thermal equilibrium by numerically studying the case of a 
free massive scalar field in de Sitter and radiation dominated cosmologies. It 
is found that a suitably defined non-equilibrium entropy per mode increases 
linearly with comoving time in a de Sitter cosmology, whereas it is not a 
monotonically increasing function in the radiation dominated case. 
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I. INTRODUCTION AND MOTIVATION 



Since its creation more than a decade ago, the inflationary universe scenario has 
become an integral part of the standard model of cosmology. However, to a great extent, 
this scenario is incomplete. One of the problems is the lack of a model of inflation that 
is adequate both from the inflationary and the particle physics viewpoints. While many 
models exist that do all the things an inflationary model needs to do, such as inflating, 
ending inflation gracefully, reheating the universe and generating safe density fluctuations, 
none of these models really is part of any reasonable extension of the standard model of 
particle physics. 

What we address in this work, however, is a more serious problem. This has to do 
with the dynamics of inflation, and more generally, the dynamics of scalar fields in an 
expanding universe. By and large, the various models of inflation make the assumption 
that the dynamics of the spatial zero mode of the (so-called) inflaton field is governed by 
some approximation to the effective potential which incorporates the effects of quantum 
fluctuations of the field. Thus the equation of motion is usually of the form: 

+ 3-0 + Kff(0) = O (1.1) 
a 

The problem here is that the effective potential is really only suited for analyzing static 
situations; it is the effective action evaluated for a field configuration that is constant in time 
0. Thus, it is inconsistent to use the effective potential in a dynamical situation. Notice 
that such inconsistency appears for any inflationary scenario (old, new, chaotic, ...). 

More generally, the standard methods of high temperature field theory are based on an 
equilibrium formalism [^,0]; there is no time evolution in such a situation. Such techniques 
preclude us from treating non-equilibrium situations such as surely exist for very weakly 
coupled theories in the early universe. 

In this work we try to rectify this situation by addressing three issues: a) obtaining the 
evolution equations for the order parameter including the quantum fluctuations, b) studying 
departures from thermal equilibrium if the initial state is specified as a thermal ensem- 
ble, c) understanding the renormalization aspects and the validity of the high temperature 
expansion. 

Our ultimate goal is to study the dynamics of phase transitions in the early universe, 
in particular, the formation and evolution of correlated domains and symmetry breaking 
in an expanding universe. From some of our previous studies on the dynamics of phase 
transitions in Minkowski space, we have learned that the familiar picture of "rolling" is 
drastically modified when the fluctuations are taken into account. As the phase transition 
proceeds fluctuations become large and correlated regions (domains) begin to grow. This 
enhancement of the fluctuations modifies substantially the evolution equation of the order 
parameter. Thus the time dependence of the order parameter is not enough to understand 
the dynamical aspects of the phase transition; it must be studied in conjunction with that 



of the fluctuations. Vilenkin and Ford and Linde |jT6[ studied the fluctuations in a free 
scalar field theory in de Sitter space and Guth and Pi [T^ studied the growth of fluctuations 
by approximating a broken symmetry situation with an inverted parabolic potential in de 
Sitter space. However, we are not aware of any previous attempt to incorporate the growth 
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of fluctuations (arising from the non-linearities) in the dynamics of the order parameter 
during cosmological phase transitions. 

Our approach is to use the functional Schrodinger formulation, wherein we specify the 
initial wavefunctional \E'[$(7);t] (or more generally a density matrix p[$(7), $(7); t]), and 
then use the Schrodinger equation to evolve this state in time. We can then use this state to 
compute all of the expectation values required in the construction of the effective equations 
of motion for the order parameter of the theory, as well as that for the fluctuations. The 
Schrodinger approach has already been used in the literature at zero temperature [1^ and 
to study non-equilibrium aspects of field theories [rT|. 



One advantage of this approach is that it is truly a dynamical one; we set up initial 
conditions at some time to by specifying the initial state and then we follow the evolution 
of the order parameter 0(t) = ($(x)) and of the fluctuations as this state evolves in time. 
Another advantage is that it allows for departures from equilibrium. Thus, issues concerning 
the restoration of symmetries in the early universe can be addressed in a much more general 
setting. 

There have been several attempts |r2|-p^] to obtain the evolution equations in expanding 
cosmologies. Our motivations, goals and many technical aspects differ substantially from 
those of previous treatments. In particular, we not only obtain the evolution equations for 
the order parameter to one-loop approximation, but we also find them in a non-perturbative 
self-consistent Hartree approximation. Within these approximation schemes, we obtain the 
evolution of the fluctuations (quantum and thermal), departures from equilibrium, study in 
detail the subtle aspects of renormalization and the validity of a high-temperature expansion. 
Our analysis applies quite generally to any arbitrary spatially flat FRW cosmology. We also 
determine the time evolution of the initial (Boltzmann) distribution functions, relate the 
time evolution to "squeezed states" and perform a numerical integration in the case of free 
fields for de Sitter and radiation dominated cosmologies. We expect to provide a numerical 
analysis of the evolution of the order parameter and the dynamics of phase transitions for 
interacting fields in a forthcoming article. 

In the next section, we set up the formalism for determining the dynamics of the order 
parameter 0(t). This involves constructing the order h equations of motion for (j){t) incorpo- 
rating quantum fluctuations, and then constructing an ansatz for the time evolved density 
matrix we need to use to evaluate the various expectation values in the problem. We then 
consider a self-consistent (Hartree) approximation to the equations of motion (sec. 3) and 
deal with the issue of renormalization of these equations (sec. 4). 

The initial state we pick for the field $(x, t) is that corresponding to a thermal density 
matrix centered at 0(t). It is then useful to try to understand the high temperature limit 
of our calculations. We are able to compute both the leading and subleading terms in the 
high T expansion of (0^(t)). From this we show that the high T expansion cannot be valid 
for all time, but breaks down in the large time limit (section 5). 

In section 6 we compute the time evolution of the Boltzmann distribution function (ini- 
tially specified as thermal) as a function of time, and find that to one-loop order and in 
the Hartree approximation, the density matrix describes a "squeezed" state. Section 7 pro- 
vides a numerical analysis of the departure from equilibrium in the simpler case of a free 
massive scalar field in de Sitter and radiation dominated cosmologies. We point out that a 
coarse-grained entropy used in the literature is not a monotonically increasing function of 
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time in the radiation dominated case. Section 8 contains our conclusions. There are two 
appendices: the first contains some technical results that are necessary in order compute the 
time evolved density matrix. The second appendix treats some of the results of the paper 
in conformal instead of comoving time. 



II. EVOLUTION EQUATIONS 

We start by setting up the Schrodinger formalism for spatially flat FRW cosmologies. 
Consider a scalar field in such a cosmology where the metric is: 

ds^ = dt^ - o^{t)dx' (2.1) 

The action and Lagrangian density are given by 



S 



j d'^xC (2.2) 



(2.3) 



Vm -\[m'-r ^7^]$2(f , t) + ^^\x, t) (2.4) 

^-4^^) (2.5) 
with TZ the Ricci scalar. The canonical momentum conjugate to $ is 

n{x,t)^a^{t)^{x,t) (2.6) 

and the Hamiltonian becomes 

= I {2^ + ^(V$)^ + a%t)V{^)^ (2.7) 

In the Schrodinger representation (at an arbitrary fixed time to), the canonical momentum 
is represented as 

r 

n(f) = -ih- 



(5$(f) 

Wave functionals obey the time dependent functional Schrodinger equation 

ih^^^^H^^^t] (2.8) 

Since wc shall eventually consider a "thermal ensemble" it is convenient to work with 
a functional density matrix p with matrix elements in the Schrodinger representation 
p[$(7), $(T);t]. We will assume that the density matrix obeys the functional Liouville equa- 
tion 
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th^=[Hit),p] (2.9) 

whose formal solution is 

p{t) = U{t,to)p{to)U-\t,to) 

where U{t, t^) is the time evolution operator, and p{to) the density matrix at the arbitrary 
initial time to- 

The diagonal density matrix elements p[^,^;t] are interpreted as a probability density 
in functional space. Since we are considering an homogeneous and isotropic background, 
the functional density matrix may be assumed to be translationally invariant. Normalizing 
the density matrix such that Trp = 1, the "order parameter" is defined as 

Ht) = ^J rf'x($(f,t)) = ^J d^xTrp{t)^x) = ^J d''xTrp{to)U-\t,toMx,to)U{t,to) 

(2.10) 

where Q is the comoving volume, and the scale factors cancel between the numerator (in 
the integral) and the denominator. Note that we have used the fact that the field operator 
does not evolve in time in this picture. The evolution equations for the order parameter are 
the following 



dt nJ ^''6^{x 

It is now convenient to write the field in the Schrodinger picture as 

^{x) = (f){t)+7]{x,t) (2.13) 
{V{x,t)) = (2.14) 

Expanding the right hand side of ( p.l2|) we find the effective equation of motion for the 
order parameter: 

d'^(f)(t) ait) debit) ,, , V"'i<\)it)) f . , 

where primes stand for derivatives with respect to 0. To leading order in the loop expansion 
we need that ijf{x,t)) = 0{h). This will be guaranteed to this order if the density matrix 
is assumed to be Gaussian with a covariance (width) 0{l/h). If ( |2.13D is introduced in the 
Hamiltonian, we arrive at: 

h'^ a{t) ^2 3 / 1, 



(2.16) 
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Keeping only the terms quadratic in rj in (p.l6| ) gives the first order term in the loop 
expansion. 

It is convenient to introduce the discrete Fourier transform of the fields in the comoving 
frame as 

V{x,t) = ^Y.vdt)e~'''' (2-17) 

In this representation, the quadratic approximation to the Hamiltonian (|2.16|) becomes the 
Hamiltonian of a collection of independent harmonic oscillators for each mode k 

= na%t)V{(f){t)) + 

^ E {-^^ + 2a^t)nimv^, + .m,,,_,] (2.18) 
Vim) = ^'(0W)v^%,o 

ujl{t) = a{t)k^ + a%t)V"{(f){t)) (2.19) 

We propose the following Gaussian ansatz for the functional density matrix elements in 
the Schrodinger representation 



k 



^vk{t)v-k{t) + ^mv-kit) + ^vk{t)fj-k{t) 



+ ^nkit) (rim - V-kit))] (2.20) 

r/fc(t) = $fe-(/.(t)\/?]%o (2.21) 

Vkit) = ^k - Ht)^h,o (2-22) 

where 0(t) = ($(a;)) and iTkit) is the Fourier transform of (n(x)). This form of the density 
matrix is dictated by the hermiticity condition $, t] = p*[$, as a result of this, 
Bk{t) is real. The kernel -Bfc(t) determines the amount of "mixing" in the density matrix, 
since if -Ba: = 0, the density matrix corresponds to a pure state because it is a wave functional 
times its complex conjugate. 

In order to solve for the time evolution of the density matrix (|2.9|) we need to specify 



the density matrix at some initial time to- It is at this point that we have to assume 
some physically motivated initial condition. We believe that this is a subtle point that 
has not received proper consideration in the literature. A system in thermal equilibrium 
has time-independent ensemble averages (as the evolution Hamiltonian commutes with the 
density matrix) and there is no memory of any initial state. However, in a time dependent 
background, the density matrix will evolve in time, departing from the equilibrium state 
and correlation functions or expectation values may depend on details of the initial state. 

We will assume that at early times the initial density matrix is thermal for the modes 
that diagonalize the Hamiltonian at to (we call these the adiabatic modes). The effective 
temperature for these modes is ksTo = l//9o- It is only in this initial state that the notion of 
"temperature" is meaningful. As the system departs from equilibrium one cannot define a 
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thermodynamic temperature. Thus in this case the "temperature" refers to the temperature 
defined in the initial state. 

The initial values of the order parameter and average canonical momentum are 0(to) = 0o 
and 71 (to) = TCo respectively. Defining the adiabatic frequencies as: 



(2.23) 



we find that the initial values of the time dependent parameters in the density matrix (p.20|) 
are 



Ak{to) 
Bkito) 



AliQ = Wk{to)a\to) coth [(3onWk{to)\ 
Wk{to)a\to) 
sinh [PonWkito)] 



Wk{to)a\to) 
4>o ; rr(to) = Ho 



tanh 



(3ohWk{tc 



(2.24) 
(2.25) 

(2.26) 
(2.27) 



Since time evolution is 
= the density matrix 



The initial density matrix is normalized such that Trp{to) = 1. 
unitary such a normalization will be constant in time. For To 
describes a pure state since = 0. 

As an example, consider the case of de Sitter space. The scale factor is given by a{t) = 
ttoC^^ and for To ^ 0, to — — oo we recognize the ground state wave-functional for the 
Bunch-Davies vacuum |10.18|. For To 7^ this initial density matrix corresponds to a 



thermal ensemble of Bunch-Davies modes. Certainly this choice is somewhat arbitrary but 
it physically describes the situation in which at very early times the adiabatic modes are in 
local thermodynamic equilibrium. Whether or not this situation actually obtains for a given 
system has to be checked explicitly. In the cosmological setting, the nature of the initial 
condition will necesarily have to result from a deeper understanding of the relationship 
between particle physics, gravitation and statistical mechanics at very large energy scales. 

Although we will continue henceforth to use this thermal initial state, it should be 
emphasized that our formalism is quite general and can be applied to any initial state. 

In the Schrodinger picture, the Liouville equation (p.9|) becomes 



ih 



dp[^,^,t] 
dt 



E 



6^ 



6^ 



2a^{t) \5rik6r]^k 6rikSr]-k, 



W{t)V'^k{m) iVk - fjk) + l^lit) {r],^r]^^ - fjkV-k) I pI'^, ^, t] 



(2.28) 



Since the modes do not mix in this approximation to the Hamiltonian, the equations for the 
kernels in the density matrix are obtained by comparing the powers of rj on both sides of 
the above equation. We obtain the following equations for the coefficients: 
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Mk 2a%t) 



(2.29) 



7 



iAk 
iBk 



Al-Bl 



aitf 



IT 



a3(t) 



(2.30) 

(2.31) 

(2.32) 
(2.33) 



The last two equations are identified with the classical equations of motion for the order 
parameter ( |2.15| ). The equation for Bk{t) reflects the fact that a pure state 5^ = remains 
pure under time evolution. 

Writing Ak in terms of its real and imaginary components Ak{t) = Ankit) + iAik{t) (and 
because Bk is real) we find that 



ARk{t) AnkiQ 
and that the time evolution is unitary (as it should be), that is 



(2.34) 



(Ankit) + Bkit)) 



constant 



(2.35) 



The initial conditions ( |2.24| , |2y25| ) and the invariance of the ratio ( |2.34| ) suggest that the 
solution for the real part of A and for B may be obtained by introducing a complex function 
Ak{t) = ARk{t) + iAik{t) 



Aukit) = ARkit) tanh [PohWkito)] = -Bk{t) sinh [l3JiWk{to)\ 
ARk{to) = Wk{to) a' (to) 
Aik{t) = Ajk{t) 

In this form, the real and imaginary parts of A satisfy the equations 

2 



ARkit) 

-Aikit) 



a3(t) 
1 



ARkit) Aikit) 

^Alkit) - A'jkit) - ujlit)a^t) 



(2.36) 
(2.37) 
(2.38) 



(2.39) 
(2.40) 



These two equations may be combined in one complex equation for the combination Akit) 
ARkit) + i Aikit) that obeys the Riccati-type equation 



iAkit) 



a%t) 



Ak^it)-ujlit)a\t) 



with the initial conditions: 



ARkito) = Wkito)a^ito) 
Aikito) = 



(2.41) 



(2.42) 
(2.43) 
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The Riccati equation ( p.41|) becomes a more amenable differential equation by the change 
of variables 



Ak{t) = -%ci\t) 



(2.44) 



The solution to the Ricatti equation with the above initial conditions is detailed in 
appendix A. We find that it is convenient to introduce two real mode functions (for each 
wavevector k) and write 

Uklit) + lUk2{t) 



(2.45) 



These mode functions obey a Schrodinger-like equation. The initial conditions on ifk{to) are 
(see appendix A) 

VkiQ = 



'Pk{t) 



\to 



ARk(t) is given by equation ( |A14| ) in appendix A, so that: 

1 



ARk{t) 

Bk{t) 



Vk{t) 



■ coth [l3ohWk{to) 



1 



^k{t) 



sinh [(3onWk{to)] 



The equal time two-point function for the fluctuation becomes 

h h 
r 



{Tlk{t)r]-k{t)) 



ipkit) \^ coth [l3ohWk{to)/2] 



2{ARk{t) + Bk{t)) 

The one-loop equation of motion for the order parameter thus becomes 



+ 3^0 + n0) + V^"'(0)^ 
a 2 



coth[P,hWk{to)/2] = 



(2.46) 



(2.47) 



(2.48) 



(2.49) 



(27r)3 2 

with the function (pk{t) defined in appendix A by (A7, A5) in which, to this order in h, 
only the classical solution <f)ci{t) enters. A consistent numerical solution of these equations 
to 0{h) would involve splitting (f){t) = (f)ci{t) + ^4'{i){t) and keeping only the 0{h) terms in 
the evolution equation. This will result in two simultaneous equations, one for the classical 
evolution of the order parameter and another for 0(i)(t). 

This equation of motion is clearly very different from the one obtained by using the 
effective potential. It may be easily seen (by writing the effective action as the classical 
action plus the logarithm of the determinant of the quadratic fluctuation operator) that this 
is the equation of motion obtained by the variation of the one-loop effective action. 

The static effective potential is clearly not the appropriate quantity to use to describe 
scalar field dynamics in an expanding universe. Although there may be some time regime 
in which the time evolution is slow and fluctuations rather small, this will certainly not be 
the case at the onset of a phase transition. As the phase transition takes place, fluctuations 
become dominant and grow in time signaling the onset of long range correlations @J^. 
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III. HARTREE EQUATIONS 



Motivated by our previous studies in Minkowski space ||,^ which showed that the growth 
of correlation and enhancement of fluctuations during a phase transition may not be de- 
scribed perturbatively, we now proceed to obtaining the equations of motion in a Hartree 
approximation. This approximation is non-perturbative in the sense that it sums up in- 
finitely many diagrams of the cactus-type |0. The Hartree approximation becomes exact 
in the N —>■ oo limit of an 0{N) vector theory. Although its validity is not warranted 
in the present case, it at least provides a consistent non-perturbative framework in which 
correlations and fluctuations can be studied. It is conceivable that this approximation could 
be implemented beyond the lowest (cactus) order in a consistent fashion. 

Vilenkin has previously studied a simplified version of the Hartree approximation in 



which, however, the mode functions that enter in the propagators did not incorporate the 
fluctuations in a self-consistent manner. 

The Hartree self-consistent approximation is implemented as follows. We decompose the 



field as in ( p.l3|) , using a potential as in (|2.4|) . We find that the Hamiltonian becomes 

+ T^W + ^Ar^^H (3.1) 



3! " ' 4! 

The Hartree approximation is obtained by assuming the factorization 

r]%x,t)^3{r]\x,t))r]ix,t) (3.2) 
r]\x, t) 6(r/2(f , t))ri'(x, t) - 3{r]\x, t)f (3.3) 

where (■ ■ ■) is the average using the time evolved density matrix. This average will be deter- 
mined self-consistently (see below). Translational invariance shows that (?7^(x, t)) can only 
be a function of time. This approximation makes the Hamiltonian quadratic at the expense 
of a self-consistent condition. In the time independent (Minkowski) case this approximation 
sums up all the "daisy" (or "cactus" ) diagrams and leads to the self-consistent gap equation 
0. In this approximation the Hamiltonian becomes 

H = na^{t)V{(f)) + 

V(0) = V(0) - hivr (3.5) 

V«(0) = n0) + ^A0(r/2) (3.6) 

V(2)(0) = n0) + ^A(r/2) (3.7) 

We can now introduce the Fourier transform of the field as in ( ^.17| ).The Hamiltonian 
will have the same form as ( p.l8| ) but the time dependent frequencies ( |2.19|) and linear term 
in 7] become 
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cul{t)=a{t)k' + a'{t)V^'\m 



(3.8) 
(3.9) 



The ansatz for the Gaussian density matrix is the same as before (|2.2CI|) , as are the evolution 
equations for the coefficients Ak{t), Bk{t), Mkit). However, the frequencies are now given 
by ( p.8|) . The c/asszca/ equations of motion (|2.32| , |2.33D now become 



-71 = V^^\<l){t))a\t) 



a3(t) 



(3.10) 
(3.11) 



The equations for the coefficients Ak{t), Bk{t), Mkit) are again solved in terms of the 
mode functions given in appendix A but with V"{(f)ci{t)) now replaced by V^'^\(j){t)) and the 
following replacement of the adiabatic frequencies: 



a{to) 



(3.12) 
(3.13) 



The mode functions of the appendix obey the differential equations and initial conditions 
with this replacement and 

I Mt) I'^l Vkit) I' • 

The initial conditions at to are now given in terms of these new adiabatic frequencies (see 
appendix A). The equal time two-point function thus becomes 



{v\x,t)) 



(2^ 



<^f(t) \^coth[(3onWk{to)/2], 



(3.14) 



which leads to the following set of self-consistent time dependent Hartree equations: 



a • 



> + 3-0 + V'(0) + A0- 
a 2 

rf2 a{t) d P 



h f d'k I ypf (t) 



(27r)3 2 

v"m)+ 



coth[/3„;i>Vfc(t«)/2] = 



^2 



n f d'k I v^it) 



(2vr) 



■coth[(3ohWk{to)/2] 



^^{t) = 



V>kit) \t= ^ 



a\to)m{tc 

m{to) 

\ a' (to) 



(3.15) 

(3.16) 

(3.17) 
(3.18) 

(3.19) 
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IV. RENORMALIZATION 



In either the one-loop approximation or the Hartree self-consistent approximation, the 
renormalization aspects are contained in the momentum integrals of the mode functions. 
Because the Bose-Einstein distribution functions are exponentially suppressed at large mo- 
menta, the finite temperature contribution will be convergent and we need only address 
the zero temperature contribution. The study of renormalization is more conveniently per- 
formed in terms of the following mode function satisfying a Schrodinger-like equation (see 
Appendix A): 



wf(t) = ai(tvf(t)yn4(^ 



3d(to) 
2a{to) 



(4.1) 
(4.2) 

(4.3) 



The one-loop case may be obtained by making the replacement yVk{4>) Wk{4>)- 

We will now analyze the renormalization aspects for the Hartree approximation; the 
one-loop case may be obtained easily from this more general case. We need to understand 
the divergences in the integral 



(27r)3 2a3(t)Wfe(to) 



(4.4) 



The divergences in this integral will be determined from the large-k behavior of the mode 
function that is a solution to the differential equation obtained from ( p.l7| ) with the initial 
conditions ( [4.2| , [4.3|) . The large-k behavior of this function may be obtained in a WKB 
approximation by introducing the function Vk{t), 



Vk{t) = exp 
Vkito) = 1 



to 



R{t')dt' 



satisfying the differential equation 



d 3 / d Id 
d^~ 2[a^2^ 



Vk{t) = 



(4.5) 
(4.6) 



(4.7) 



with irf^x^t)) being the self-consistent integral in the Hartree equation ( |3.17| ). The one- 
loop approximation is obtained by setting {rf{x,t)) = in the above equation. The mode 
function is obtained as a linear combination of the function Vk{t) and its complex 
conjugate; the coefficients are to be determined from the initial conditions (see below). The 
function R{t) obeys a Riccati equation 



■ „9 3 / d 1 d . 
R + R^ --{- + -—] + 



2 \a ' 2a2 

We propose a WKB solution to this equation of the form 



(4.8) 
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R 



a{t) °' ' k 

and find the time dependent coefficients by comparing powers of k. Tliis yields: 

d{t) 



Ro{t) 
Ri{t) 



2a{t) 
a{t) 
~Y 
1 d 



Finally, we write 



i^^it) = l[i+i]m) + l[i-imt) 



where 7 is determined from the initial condition (f4.3|) : 



7 = 1 — ^- 



+ 



k ' 2P 
Thus in the k 00 limit we find 

KWP f 1 , 1 \d\Q 



f 1 j_ 



Introducing an upper momentum cut-off A we obtain 



(2vr) 



h /A 



2 

a2(t) 



coth [PonWkito)/2] 



h A2 



57r 



2a2(t) 



+ 



-| + \/"(0(t)) + ^(r/2(f,t)) 



finite 



(4.9) 



(4.10) 
(4.11) 
(4.12) 

(4.13) 
(4.14) 



+ Oil/k") + 



(4.15) 



(4.16) 



where we have introduced a renormalization point K, and the finite part depends on time, 
temperature and K. In the one- loop approximation, {7f{x,t)) does not appear in the log- 
arithmic divergent term in ( [4.16| ) (as it does not appear in the differential equation for the 
mode functions up to one loop). 

There are several physically important features of the divergent structure obtained above. 
First, the quadratically divergent term reflects the fact that the physical momentum cut-off is 
being red-shifted by the expansion. This term will not appear in dimensional regularization. 

Secondly, the logarithmic divergence contains a term that reflects the initial condition 
(the derivative of the expansion factor at the initial time to)- The initial condition breaks 
any remnant symmetry. For example, in de Sitter space there is still invariance under the de 
Sitter group, but this is also broken by the initial condition at an arbitrary time to. Thus this 
term is not forbidden, and its appearance does not come as a surprise. As a consequence of 
this term, we need a time dependent term in the bare mass proportional to l/a^{t). 
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We are now in a position to present the renormalization prescription within the Hartree 
approximation. In this approximation there are no interactions, since the Hamiltonian is 
quadratic. The non-hnearities are encoded in the self-consistency condition. Because of 
this, there are no counterterms with which to cancel the divergences and the differential 
equation for the mode functions ( |3.17|) must be finite. Thus the renormalization conditions 
are obtained from 



(4.17) 



where the subscripts B, R refer to bare and renormalized quantities respectively and {ri'^)B 
is read from ( [4.16|) 



A2 h , 

+ — ^In 



Sir^a^it) 87r2 



a\to) 
a^{t) 



-^+ml + ^<P'{t)+^nn+^^'-' 



n 

6" 



R 



Using the renormalization conditions ( [4.17| ) we obtain 



m|(t) + 
\b = — 



Xsh Xb^ 
167r2a2(t) ^ 16^ 

Xr 



In — 

K 



a^t) 



1 



16 

Xsh 



R 



167r2 
Ir + J 



ln(A 
K 



IR 



1 + ^ln 

167r2 



+ finite 
(4.18) 



(4.19) 
(4.20) 

(4.21) 
(4.22) 



where 



IR 



(2 



J 



e{k - K) 

4P 



+ 



n 

'~6 



_ a \to] 
a\t) 



2ka^{t) 

■+< + ^^\t)+iRT^+'-^{v')R 



Xr 
2 



(27r)3exp/3„;iWfe(to)-l 



(4.23) 
(4.24) 



The conformal coupling ^ = 1/6 is a fixed point under renormalization [jT8[. In dimen- 
sional regularization the terms involving are absent and In A is replaced by a simple pole 
at the physical dimension. Even in such a regularization scheme, however, a time dependent 
bare mass is needed. The presence of this new renormalization allows us to introduce a new 
renormalized mass term of the form 



a?{t) 



This counterterm may be interpreted as a squared mass red-shifted by the expansion of the 
universe. However, we shall set E = for simplicity. 
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It is clear that there is no wave function renormahzation. This is a consequence of the 
approximations invoked. There is, in fact, no wavefunction renormahzation in either the 
one-loop or the Hartree approximation for a scalar field theory in three spatial dimensions. 

Notice that there is a weak cut-off dependence on the effective equation of motion for 
the order parameter. 

For fixed X^, as the cutoff A ^ oo 



In addition, 



A; 



(Any 



(4.25) 



mlit) 



1 



In (A 



+ 



In A 



(4.26) 
(4.27) 



This approach to 0~ of the bare coupling as the cutoff is removed translates into an 
instability in the bare theory. This is a consequence of the fact that the N-component $^ 
theory for — oo is asymptotically free (see ref. pO|), which is not relieved in curved 
space-time. Clearly this theory is sensible only as a low-energy cut-off effective theory, and 
it is in this restricted sense that we will ignore the weak cut-off dependence and neglect the 



term proportional to the bare coupling in ( [4.28 



The renormalized self-consistent Hartree equations thus become after letting A = oo: 



R 



A 



: 



(4.28) 
(4.29) 



where (?7^)_r is given by equations ( |4.23| , |4.24| ). 

For completeness we quote the renormalized equation of motion for the order parameter 
to one-loop 



a ■ 



+ 3-0 + m^0 + ei^7^0 + ^0=^ + ^{v') 



R 







d{t) d k'^ 
d^^ ^Jt^'^ 



+ m^ + ei^7^+^02(^) 



ipk{t) = 



(4.30) 

(4.31) 
(4.32) 



The one-loop coupling constant renormahzation differs from that in the Hartree approxi- 
mation by a factor of three. This is a consequence of the fact that the Hartree approximation 
is equivalent to a large N approximation and only sums up the s-channel bubbles. 
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V. HIGH TEMPERATURE LIMIT 



One of the payoffs of understanding the large-A; behavior of the mode functions (as 
obtained in the previous section via the WKB method) is that it permits the evaluation of 
the high temperature hmit. We shall perform our analysis of the high temperature expansion 
for the Hartree approximation. The one-loop case may be read off from these results. 

The finite temperature contribution is determined by the integral 



J = n 



d'k I (t) 



(5.1) 



For large temperature, only momenta k > To contribute. Thus the leading contribution is 
determined by the first term of the function R{t) (eq. (|4.9| )) of the previous section. We find 



J 



12h 



a{t) 



[1 + 0{l/To) 



(5.2) 



Thus we see that the leading high temperature behavior reflects the physical red-shift in 
the cosmological background and it results in an effective time dependent temperature 



Teffit) 



a(to) 
a{t) 



To leading order, the expression obtained for the time dependent effective temperature 
corresponds to what would be obtained for an adiabatic (isentropic) expansion for blackbody- 
type radiation consisting of massless relativistic particles evolving in the cosmological back- 
ground. 

This behavior only appears at leading order in the high temperature expansion. There 
are sub leading terms that must be taken into account. These can be calculated within the 
high temperature expansion and we do this below. To avoid cluttering of notation, we will 
set fc^ = = 1 in what follows. 

We define 



A 



(5.3) 



and we will assume that m'^{To) <^ T^. Since we are interested in the description of a phase 
transition, we will write 



24 ' 



Thus, to leading order in To 



T^ > 



(5.4) 



m'{To) 



24 



{T^ - T^) 



Our high temperature expansion will assume fixed m(To) and m(To)/To <^ 1. 
It becomes convenient to define the variable 



(5.5) 
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X 



Recall from our WKB analysis that the leading behavior for A; — > oo is (see equation ( |4.15| ) 

I ^Ut) ? 1 



2a^{t)k 



adding and subtracting this leading term in the integral J and performing the above change 
of variables, we have 



J = J1 + J2 
a{to) 



Ji 
J2 



ait)) 



T 



TT 



27r2 

\rp2 

2 ° 



a{to) 
a{t) 

a{to) 



^ xdx 
xdx 



a\t)Jx^T^-m^{To) 



m{To) _ Y 

To 



a{t) 



m{To) ^ m\To) 



12 27r2To 87r2T2 



O 



' m\To) 
, T3 



a{t) 
2a{to] 



(5.6) 



+ 



(5.7) 



We now must study the high temperature expansion of Ji. We will restrict ourselves 
to the determination of the linear and logarithmic dependence on Tq. For this purpose, it 
becomes convenient to introduce yet another change of variables 

m{To) 
X = — = — z 



and use the fact that in the limit To ^ ?Ti(To 

z 



e To ' 



%[l-^z + 



This yields the following linear and logarithmic terms in To'. 



J 



lUn 



a{t) 



Tom{To) 



dz \a^{t)m{To)V^^ 



a{t) 



(5.8) 



Note that the above integral is finite. 

The logarithmic contribution is obtained by keeping the 0{l/k^) in the large momentum 
expansion of | ip^{t) P given by equation ( [4.15| ) (in terms of the new variable z). We obtain 
after some straightforward algebra: 



Ji 



In 



llog 



m{To) 
To 



n d^{to) 



a2(t) 



+ 



a{to) 



a{t) 



m\T,) + 



24 



24 



(5.9) 



That is, in the limit >> m(To), Ji = Jm^ + Juog + 0((To)° 
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Comparing the 0{T^, To, In To) contributions it becomes clear that they have very differ- 
ent time dependences through the scale factor a{t). Thus the high temperature expansion as 
presented will not remain accurate at large times since the term quadratic in To may become 
of the same order or smaller than the linear or logarithmic terms. The high temperature 
expansion and the long time limit are thus not interchangeable, and any high temperature 
expansion is thus bound to be valid only within some time regime that depends on the initial 
value of the temperature and the initial conditions. 

As an illustration of this observation, we calculate Jmn explicitly in the case of de Sitter 
space. We need to obtain | (p^(t) p in order to evaluate the integral in ([5.81). Inserting 
the term proportional to in the Hartree equations, we find that | f^i^) P obeys the 
differential equation 



^+m\To 



+ 



M2 

24 



-2Ht 



Ml 

24 



The solution of this equation is given by: 



C2HfhBke 



'3/2 



\Ht 



(5.10) 



(5.11) 



Bk 



1 



+ m[^lo}+ 24 



(1 2) 

where H^j^ are the Hankel functions and we have assumed 



m\To) 



Ml 

24 



< H. 



The coefficients Ci, C2 are determined by the initial conditions on fki^) described above. 
We finally obtain 



llin 



m(To)To 



°° dz -J - 1 



ArT2 
24 



^ + m'^{T)z 



(5.12) 



This term is time independent, finite and positive. This example clearly illustrates the fact 
that different powers of To enter in the expansion with different functions of time and that 
the high temperature expansion is non-uniform as a function of time. 



VI. EVOLUTION OF THE INITIAL DISTRIBUTION 



The initial density matrix at t = to was assumed to be thermal for the adiabatic modes. 
This corresponds to a Boltzmann distribution for the uncoupled harmonic oscillators for the 
adiabatic modes of momentum k, and frequencies Wfc(to) (in the Hartree approximation; as 
usual, the one-loop result can be found by replacing this with Wk{to))- That is 



Ho = J2 ^^k{to) al{to)ak{to) + 



(6.1) 
(6.2) 
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The creation and annihilation operators define the initial occupation number of the adiabatic 
modes: 



(Mto)) ^ 



(6.3) 
(6.4) 



where the expectation value in (6^) is in the initial density matrix at time to- 

In a time dependent gravitational background the concept of particle is ill-defined. How- 
ever, by postulating an equilibrium initial density matrix of the above form, a preferred 
"pointer" basis is singled out at the initial time. It is this basis that provides a natural 
definition of particles at the initial time and we can use it to ask: how does the expectation 
value of this number operator evolve in time? 

At any time t, this expectation value is given by 



{Nk){t) 



Tral{to)ak{to)p{t) 
Trp{to) 



(6.5) 



This quantity gives information on how the original Boltzmann distribution function for 
the adiabatic modes evolves with time. The k = mode will receive a contribution from 
the order parameter, but since the number of particles is not conserved (no charge) there 
is no bose condensation and the k = mode will give a negligible contribution to the total 
number of particles. Thus we only concentrate on the k ^ modes. 

The expectation value ( |6.5|) may be easily computed by writing the creation and anni- 
hilation operators in terms of r^^, 11^ = 5/5rj_k in the Schrodinger picture at to- The result 
of doing this is: 



al{to) 

Oikito) 



1 



2h 
1 



1 



a\to)Wk{to) % 
1 5 



+ \j a\to)Wk{to)7]--k 



a^{to)mito)vk 



After some straightforward algebra we find 

{Nk){t) + ^ = (2 I Mt,to) I' -l) ((iV,)(t„) + i 

,2 1 I Vfit) P L , a^t) I 0f (t) p 



(6.6) 
(6.7) 



^k(t, to) 



4 I ^^(t) P o 



1 + 



a^ito)Wlito) I v^it) 



+ 



(6.8) 
(6.9) 



where we have made use of ( |A14| , |A15|) . This result exhibits the two contributions from 
"spontaneous" (proportional to the initial thermal occupation) and "induced" (independent 
of it). 

We now show that this result may be understood as a Bogoliubov transformation. To 
do this, consider the expansion of the field in the Heisenberg picture: 



Vkit) = ^{ak^^\t) + aiipj!it) 



(6.10) 
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where the mode functions satisfy 



a dt 



(6.11) 



together with the self-consistency relation. Then the Heisenberg field rik{t) is a solution of 
the Hartree Heisenberg equations of motion and the al , at create and destroy the Hartree- 
Fock states. Notice that in the Heisenberg picture, these creation and annihilation operators 
do not depend on time. Using the Wronskian properties of the functions ^]}{t) (see appendix 
A) we can invert and find the creation and annihilation operators in terms of rik{t) and its 
canonically conjugate momentum H_fc(t). Once we have expressed these operators in the 
Heisenberg picture in terms of the field and its canonically conjugate momentum, we can go 
to the Schrodinger picture at time to- In this picture the creation and annihilation operators 
depend on time and are given by 



(6.12) 
(6.13) 



^k{t) = ^ \ak{to)^^{t)-a\t)r^k{to)^Ut) 
al,{t) = ^ [H,(t„)¥,f (t) - a\t)r^k{to)Vk\t) 

The Schrodinger picture fields at to can be written in terms of the operators ( |6.6| , 
and we finally find the creation and destruction operators at time t to be related to those 
at time to by a Bogoliubov transformation: 



(6.14) 



If we now compute the average of the new creation and annihilation operators in the initial 
density matrix and write the mode functions (pk{t) in terms of the real functions Ui^2 as 
defined in appendix A, we recognize | JF+ kit, to) p to be the same as | J^k(t,to) ^ given by 

(PD- 

We also find that 



^ +,A:(^, to) 



^-,k{t, to) 1^— 1 



as is required for a Bogoliubov transformation. 

One way to interpret this result is that, at least within the one-loop or Hartree approxi- 
mations, time evolution corresponds to a Bogoliubov transformation. This interpretation is, 
in fact, consistent with the result that in these approximation schemes, the density matrix 
remains Gaussian with the only change being that the covariance and mixing terms change 
with time. 

Thus within the one-loop or Hartree approximation, time evolution corresponds to a 
"squeezing" of the initial state. The covariance changes with time and this corresponds to 



a Bogoliubov transformation. As argued by Grishchuk and Sidorov |2T|] the amplification of 
quantum fluctuations during inflation is a process of quantum squeezing and it corresponds 
to a Bogoliubov transformation. The properties of these "squeezed" quantum states have 
been investigated in references P^j-E^ . 
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Hu and Pavon ||25|, Hu and Kandrup p6|] and Kandrup have introduced a non- 
equilibrium, coarse-grained entropy that grows in time as a consequence of particle pro- 
duction and "parametric amplification" . This definition was generalized by Brandenberger 
et al. and Gasperini and collaborators to give a measure of the entropy of the 



gravitational field. The growth of this entropy is again a consequence of the parametric am- 
plification of fiuctuations and the "squeezing" of the quantum state under time evolution. 

These authors argue that the non-equilibrium coarse-grained entropy in the mode of 
( CO moving) wavevector k is 



(6.15) 



in the case when {Nk){t) ^ 1. Thus the growth of entropy is directly associated with 
"particle production" or in our case to the evolution of the initial Boltzmann distribution 
function. 

The coefficient of parametric amplification is related to the Bogoliubov coefficient given 
by equation (|S.9|). Thus this coefficient directly determines the time dependence of the 
non-equilibrium coarse-grained entropy. 



VII. TWO SIMPLE EXAMPLES: 

As an example of the method that allows the out-of equilibrium time evolution, we solve 
numerically the simple cases of a free massive scalar field in two relevant cosmologies: 

A. de Sitter cosmology: 

In this case the scale factor is a{t) = OoC^* with H being Bubble's constant. The impor- 
tant quantities that encode the out of equilibrium evolution are the mode functions <fk{t) 
that obey the equation ( [A^ ) with V"{(f)) = m?. With these mode functions we construct the 
real and imaginary parts of Ak{t) and the kernels of the density matrix (see eq. |2.36|) and 



the parametric amplification factor ( |6.9| ). It is convenient to rescale the differential equation 
and functions and define 



, ke m 

z = Hit-Q;q = ——;e = - 



(7.1) 



the variable q is recognized as the physical wavevector at the initial time to multiplied by 
the horizon size (or equivalently, horizon size divided by physical wavelength). At the same 
time it is convenient to rescale the mode functions (in terms of their real and imaginary 
parts) 



R,I/ 



^f(O) 



^o(O) = 



dz 



\z=0 



dz 



\z=0- 



2 I 2 

q + e 



(7.2) 
(7.3) 

(7.4) 
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The differential equation becomes 



dz 



R,Ii 







(7.5) 



We have numerically integrated these equations with the above boundary conditions choosing 
as representative parameters m = 1 Gev ; H = 10^° Gev, q = 0.01, 0.1, 1, 10, 100. In figure 
(l.a,b,c) we show ARkl{.o?^H) for q = 0.1,1.0,10 as function of 2; for e = 10~^°. Figure 
(2) shows the logarithm of Aike'^^*" /{alH) for q = 0.01,0.1,1.0,10,100 function of 
z. Whereas the real part tends to (a g-dependent) constant at large times, the imaginary 
part grows at long times as ~ ^H{t~to) ^^j^g slope on the graph is 1). Figure (3) shows the 
logarithm of the parametric amplification factor ( |6.9| ) for the above values of g as a function 
of z. The slope of the lines at long times is 2, thus the parametric amplification factor 



grows as 



,2H{t-to) 



at large times. Thus we clearly see that if the distribution function 
at some initial time to was determined by an equilibrium Boltzmann factor (in terms of 
the comoving wavelengths), this distribution function evolves in time out of equilibrium 
and grows with time approximately as ~ g2H(t-to) long times but with different rates for 
different wavevectors. 

From the numerical integration we see that in the case of de Sitter expansion, at late 
times (when the number of "particles produced" is large and the expression for the entropy 
( |6.15| ) is valid) the entropy per mode grows linearly with comoving time. 

B. Radiation dominated cosmology 



For a radiation dominated cosmology the scale factor is given by 

a{t) = ao (^Y^ 

For this convenient rescaling is in terms of the variables 

k 



s = mt ] q 



mao 

The differential equation for the mode functions, and boundary conditions thus become 

(7.6) 



+ ^ ^ + ^' + 1 
ds'^ 2s ds s 







+ 1 







g' + l 



(7.7) 
(7.8) 



ds ds 
There are several noteworthy features in the radiation dominated case as exhibited 
by figures (4.a-d) and (5.a,b). There is a strong dependence on the initial condition as 
parametrized by to {so = mto) and also a strong dependence on the initial physical wavenum- 
ber q. Perhaps the most notable feature are the oscillations in both the real and imaginary 
parts of the covariance A that translate into an oscillatory behavior in the parametric ampli- 
fication factor. As result, at long times the entropy per mode ( |6.15|) is not a monotonically 
increasing function of comoving time. 
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VIII. CONCLUSIONS 



Non-equilibrium aspects of the dynamics of scalar fields in spatially flat FRW cosmolo- 
gies were studied by means of a functional Schrodinger approach. The initial state was 
specified as a thermal density matrix at some early initial time assuming local thermody- 
namic equilibrium for the adiabatic modes at that particular time. This density matrix was 
evolved in time and the evolution equations for the order parameter (ensemble average of the 
scalar field) and the fluctuations were obtained both to one-loop and in a non-perturbative 
self-consistent Hartree approximation. The renormalization aspects were studied in detail 
and it was pointed out that the renormalization counterterms contain a dependence on the 
initial conditions through the scale factor and its derivatives at the initial time. 

The high temperature expansion was investigated and it was found that the limit of high 
temperatures and long times are not interchangeable. As a consequence of the red-shift of 
the initial temperature the coefficients of the different powers of temperature are different 
time-dependent functions. The high temperature expansion is only valid within a short time 
interval after the initial time and certainly breaks down at long times. 

The time evolution of the Boltzmann distribution functions (initially the thermal equilib- 
rium distribution functions) is obtained. It is pointed out that to one-loop order and also in 
the Hartree approximation, the time evolved density matrix describes quantum "squeezed" 
states and time evolution corresponds to a Bogoliubov transformation. 

To illustrate the departure of equilibrium, we have studied numerically the case of a 
free massive scalar field in de Sitter and radiation dominated cosmologies. It was found 
that a suitably defined coarse-grained non-equilibrium entropy (per k mode) grows linearly 
with time in the de Sitter case but it is not a monotonically increasing function of time 
in the radiation dominated case. This result may cast some doubt on the applicability of 
this definition of the non-equilibrium entropy. There still remain some (open) fundamental 
questions regarding the connection of this entropy and the thermodynamic entropy of the 
universe, in particular whether the amount of entropy produced is consistent with the current 
bounds. 

This work sets the stage for a numerical study of the dynamics of phase transitions 
in cosmology fully incorporating the non-equilibrium aspects in the evolution of the order 
parameter and which at the same time can account for the dynamics of the fluctuations 
which will necesarily become very important during the phase transition. 

We expect to report on the numerical study of the phase transition in a forthcoming 
article [|8[. 
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APPENDIX A: 



The Riccati equation ( f^.41D can be transformed into a linear differential equation by the 
change of variables 



Ak{t) = -ia\t) 



3m Vkit) 



Vk{t) 



(Al) 



We find that ^kit) obeys a simple evolution equation 



+ 3--^ + 
a dt 



(A2) 



with (j)ci(t) the solution to the classical equation of motion ( p.32| , p.33| ). In the Hartree 
approximation V'^cpciii)) should be replaced by V^'^\(f){t)) with (f)(t) the full solution of the 
self-consistent equations. 



d^^k _^^ddip^ 



dt^ 



a dt 



= 



(A3) 



The Wronskian for two arbitrary solutions to the above differential equation is 

C 

W[v9i, LP2] = ip2{t)ipi{t) - ipi{t)ip2{t) 



a3(t) 



(A4) 



with C a constant. Then writing (pk{t) = ^ik{t) + i^2k{t) with (/?i.2 real solutions we find 

C 



Aukit) 



[v\k + Vlk] 



Aik{t) = -a^{t) 



vlk + 'Ak 



(A5) 
(A6) 



It proves convenient to introduce (for all k) the real functions Ui^2{t) as 



Wk{t,) 



(A7) 



The Uak for a = 1,2 are real and satisfy the Schrodinger-like differential equation 
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(A9) 



with W[- ■ ■] the Wronskian, and C is the same constant as above. Since the choice of C 
corresponds to a choice of normahzation of these functions, we choose C = 1. The initial 
conditions ( p.42| , 2.43| ) still leave one free condition on these functions, we choose 



Uik{to) ^ 

U2k{to) = 

The boundary conditions on the mode functions Ua^k are 
Uik{to) = 1 ; U2k{to) = 



3 a{to) 



2a{to) 



U2k{to) = Wk{t, 



V"{<Pci{to)) 



(AlO) 
(All) 



(A12) 

(A13) 



and the corresponding replacement for the Hartree case. Thus the final solution to the 
Riccati equation ( |2.41| ) with the given initial conditions is 



ARk{t) 



a%t)Wk{to) 
Aik{t) = -a\t) 



Ulk [Uik — ff^lfc] + U2k U^2k — |f^2fc 



(AM) 
(A15) 



u!kit)+ul,{t) 

In terms of the original functions ipk{t) ( |Al| ) the initial conditions are simply 

^k{to) = ^ (A16) 



a'iQWkit, 



WkiQ 



(A17) 



Then the initial conditions for ipk{t) ; ^pli^) are naturally interpreted as those for negative 
and positive (adiabatic) frequency modes at the initial time to- 

The reason for introducing the functions Ua,k(t) is because these obey a simpler second 
order Schrodinger-like equation which is amenable to be studied in the asymptotic regime 
via WKB approximations (see the section on renormalization) . 



APPENDIX B: CONFORMAL TIME ANALYSIS 

It is interesting to see how some of our results can be obtained by rewriting the metric 
in terms of the conformal time defined by: 



V 



* dt' 



(Bl) 
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The first thing we should note is that the physics should not depend on what time 
coordinate is used, since the theory should be generally coordinate invariant (there are no 
gravitational anomalies in four dimensions [^). Thus, the field amplitude and the canonical 
momentum in conformal time should be related to those in comoving time via a canonical 
transformation. We now show that this is indeed the case. 

Using the conformal time version of the line element ds^ = a?{r]){drf — dx'^), the scalar 
field action can be written as: 



S[<!>]= / dr]d^xa\r]) 



((9,$)2-(V$)2)-\/($(r/,x)) 



(B2) 



where the potential term is as in the text (i.e. it could include a coupling to the curvature 
scalar). Following the standard procedure for obtaining the canonical momentum n(?7, x) to 
$(?7,a;), and to get at the Hamiltonian density yields: 

n2 a\r]) 



n 



+ 



-(V$) +a^(r7)l^(<l>). 



(B3) 



2a2(?7) ■ 2 

Here conformal time derivatives are denoted by a prime. The generator H of displacements 
in conformal time is the spatial integral of Ti above. Thus the conformal time Liouville 
equation reads: 



(B4) 



If we label the field and its conjugate momemtum in the comoving time frame (i.e. that 
of the text) as 11 respectively, the results of section 2 are: 



n(f,t) = a^{t)^{x,t) 

H= /rf3^/iI^ + ^(v|.)2 + aMl>) 



^2,3 ■ J 

The Liouville equation in comoving time can be rewritten in conformal time using the 
relation: d/dt = a{ri)~^ d / dr] . After doing this we find that eq. ( |2.9| ) becomes: 

.5p(n,^) 



a(r7)/7,/)(n,<l>) 



But 



drx 



(B6) 



(B7) 



Comparing this with eq.(^^), we see that we can make the identifications: 11 = 11, \1/ = 
Thus not only is the physics equivalent in both coordinate systems (as must have been 
the case), but the physics in the two coordinate systems related by a trivial canonical 
transformation. 

We can rewrite all of the comoving time results in terms of conformal time. Some of the 
equations take on a much simpler form in conformal time than in comoving time. This will 
be important when numerical issues are tackled, such as an analysis of the back-reaction 



problem. This work is in progress [£8 
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Figure Captions: 

Figure 1(a): Ar/ a^H vs. z for g = 0.1, m/H = 10~^°. 

Figure 1(b) : AR/alH vs. z for de Sitter for g = 1.0, m/i/ = 10"^°. 

Figure 1(c)" : AR/alH vs. z for de Sitter for q = 10, m/i/ = 10"^°. 

Figure 2 : In (^7e-^^*Va^//) vs. z for de Sitter for g = 0.01; 0.1; 1; 10; 100, 
m/H — 10~^^. The slopes of the hnes at long times is 1. 

Figure 3 : In \ \ vs. z for de Sitter for q = 0.01; 0.1; 1; 10; 100, m/H = 10-^°. 
The slope of the lines at long times is 2. 

Figure 4(a): AR/a^m vs. s for radiation dominated for g = 1 (solid line) and q = 10 
(dashed line); Sg = 1.0. 

Figure 4(b): Ai/a^m vs. s for radiation dominated for g = 1 (solid line) and q — 10 
(dashed line); So = 1. 

Figure 4(c): AR/a^m vs. s for radiation dominated for g = 1 (solid line) and g = 10 
(dashed line); Sq = 10. 

Figure 4(d): Ai/alm vs. s for radiation dominated for g = 1 (solid line) and q — 10 
(dashed line); So — 10. 

Figure 5(a): | ^ |^ vs. s for radiation dominated for g = 1 (solid line) and g = 10 
(dashed line); So = 1 

Figure 5(b): | p vs. s for radiation dominated for g = 1 (solid line) and g = 10 
(dashed hne); Sg = 10 
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